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Abstract 

Production of electron-positron pairs with invariant masses less than 300 MeV 
from thermalized hadronic matter in relativistic heavy ion collisions is calculated 
using a soft virtual photon approximation. The general theoretical framework is re- 
viewed and extended to include arbitrarily massed and charged reaction partners, 
which we apply to pions and quarks. This result, exact within the soft photon ap- 
proximation, is compared with a widely used approximate result which uses an elec- 
tromagnetic amplitude limited in validity to momentum transfers less than 4m^ q y 
This momentum-restricted method works very well for pions, constituent quarks and 
medium quarks; whereas, it fails when applied to current quarks. A field theory 
calculation is performed for the irir elastic cross section which gives excellent agree- 
ment with data. Quark-antiquark annihilation diagrams in the Born approximation 
are estimated and compared with ir° and r/ Dalitz decay contributions to the e + e~ 
invariant mass spectra. A comparison is made between the rate of production of 
zero total momentum soft dileptons obtained using resummation techniques in QCD 
perturbation theory to that which we calculate using this soft photon approximation. 
Then a Bjorken picture for the evolution is adopted allowing an integration over the 
history of the colliding nuclei. Using initial conditions likely to be found at RHIC or 
LHC we conclude even if Dalitz decays can be identified and subtracted from the ex- 
perimental data it will be difficult to distinguish quark from pion degrees of freedom 
in the low-mass e + e~ spectra. 
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I. INTRODUCTION 



Dileptons and photons (T^U provide natural tools for observing thermalized nuclear 
matter. Once produced, they hardly interact and collide infrequently, if at all, with hadrons 
in a region of hot and dense matter. Therefore, photons and lepton pairs are among the 
best carriers of information about the earliest stages of the heavy ion collision. In principle, 
they can carry information about the temperature of the primordial state. That there are 
many competing sources of photons and dileptons both before and after the hot stage of the 
collision makes reconstruction of the hottest state very difficult. The best that one can do is 
to evaluate the contribution from each source to ascertain whether or not it is dominant in 
some part of the phase space or if it has distinct properties rendering separation feasible. We 
proceed along these lines by considering four regions in the invariant mass M of the produced 
dileptons. For invariant masses M above the J/ if) peak, the spectrum is dominated by Drell- 
Yan production ||. The spectrum around M ~ 1 GeV is dominated by p, u, and 4> decays. 
Dileptons with invariant masses in the region 1.5-3.0 GeV seem to be mostly thermal ||. 
Competing sources in this mass region are DD decays H, Drell-Yan and pre-equilibrium 
emission Recent calculations p] of dilepton spectra at invariant mass 1 GeV < M < 

3 GeV show that for presently achieved energy densities (CERN, BNL) the background is 
larger than thermal contributions. This is somewhat contrary to earlier calculations 
and the reason, in our opinion, is that the earlier calculations used initial temperatures and 
energy densities that were too high. 

The last region of invariant mass (M < m p ), which we will consider in this paper, is a 
region of soft dilepton production ||. There are many sources of dileptons with masses M < 
7T° and rj Dalitz decays, tt + tt~ annihilation (for M e+e > 2m n ), pion scattering with 



m 



p- 



virtual bremsstrahlung and if a quark-gluon plasma (QGP) is formed there will be quark- 
antiquark annihilation and quark-quark(antiquark) scattering with virtual bremsstrahlung. 
Recent calculations ||TQ|JTT[1 have shown that in the mass region M ~ Im^ the contribution 
from rj Dalitz decay is two or thee orders of magnitude above qq annihilation spectra. But 
these calculations of the qq annihilation were done in the Born approximation. However, it 
has been shown [I2| in the case of zero p± dileptons that for small mass pairs the perturbative 
corrections can be several order of magnitude larger than the Born term. So, we consider 
this qq annihilation spectra as a lower limit only. 

The contribution of virtual bremsstrahlung from pion scattering is the same order as 
rj Dalitz decay in the mass region M ~ 2m n . This conclusion was reached by including 
only the processes tc^tt — > 7r ± 7r°7* [[K]]. There are other pion-pion reactions, such as 
7r ± 7T :F — > 7T ± 7r =F 7*, which will contribute to the overall e + e~ spectrum. We will discuss such 
processes in this paper. Typically, a momentum transfer restriction (\t\ < 4m 2 ) is the price 
to pay for a simple but approximate expression describing the electromagnetic factor in the 
cross section. Such expressions can be quite accurate when the masses are large, but for 
small massed particles such as current quarks, the approximation may not be so good. In 
this paper we relax this restriction on \t\ by considering scattering at any angles or energies 
for the charged particles. We will investigate the influence of such relaxation on the e + e~ 
production rates and total yields. 

If we consider process of qq annihilation in a QGP, it is also necessary to include quark 
(antiquark) and/or gluon scattering with virtual bremsstrahlung and their contribution to 
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pair production. Their consideration is not possible without knowledge of the differential 
cross sections for quark-quark and quark-gluon elastic scattering. Using lattice gauge theory- 
data on screening, we shall estimate the current (m q ~ 5 MeV), medium (m q = [2na s /$\ l l 2 T) 
and constituent (m q ~ 300 MeV) quark contribution to the electron-positron mass spectrum. 

The paper is organized in the following way. In Sect. |I| we present the picture of evolution 
for the colliding nuclear system. In Sect. |T| we discuss the soft photon approximation for 
general quasi-elastic scattering with virtual bremsstrahlung and calculate the rate of e + e~ 
production. The two contributing elements are the amplitude for elastic scattering and the 
polarization dotted into the electromagnetic current. We derive the exact electromagnetic 
factor and compare it with a popular approximate result, limited in validity to \t\ < 4m 2 . 
In Sect. |V| we apply this formalism to thermalized pions at temperature T. Again, the 
exact and approximate results are presented. Then in Sect. |V| we use our formulae to 
calculate pair production through virtual bremsstrahlung in quark (antiquark)-gluon and 
quark-quark (antiquark) scattering. A comparison is made in Sect. [VI] among rates of zero 
(total) momentum pairs calculated using this soft photon approximation, QCD perturbation 
theory and quark-antiquark annihilation with a modified dispersion relation. In order to 
obtain estimates for the total yields from these competing pion and quark processes, we 
integrate the rates over the time or temperature evolution of the colliding nuclei within a 
Bjorken picture |H| and present the results in Sect. |VTi| . Also in Sect. |V11| we use isentropic 
evolution approximations to estimate tt°, rj Dalitz decays and qq annihilation in the Born 
approximation. Finally, in Sect. [VIII| we write a brief summary and present concluding 
remarks. 



II. EVOLUTION OF THE COLLIDING NUCLEAR SYSTEM 

We assume formation of thermalized hadronic matter in relativistic heavy ion collisions. 
This approach has recently been popularized especially with respect to QGP formation. 
There is evidence that in such collisions complete equilibration occurs within 1-2 fm/c ||14| . 
The idea of thermalization is attractive because it enables us to use powerful methods from 
equilibrium thermodynamics. If the initial temperature Tj > T c ~ 200 MeV, then we expect 
hadrons to dissolve into freely propagating quarks, antiquarks and gluons after having gone 
through a deconfinement phase transition. It is well known [[| that chiral invariance is 
restored near the phase transition temperature, whereupon constituent quarks (m q ~ 300 
MeV) become current quarks with masses m q ~ 5 MeV. On the other hand, in a QGP the 
fermion mass induced by many-body effects is || m q = [27ro; s (T)/3] 1 / 2 T. If we adopt the 



lattice gauge theory renormalization group result [15 



" s(T) (33-2JV»lii(aT/T c )' (2-1) 

with a ~ 8, then for relevant temperatures T ~ 200-300 MeV the medium quark masses 
m q ~ 100-150 MeV. So, if we want to calculate the bremsstrahlung contribution from the 
quark phase, it is necessary to take into account the processes of current and medium quark 
scattering at temperatures T > T c as well as constituent quark scattering at T = T c . During 
a time evolution the nuclear matter cools. If the QGP-hadron phase transition is first 
order at T = T c , the mixed phase may be realized. In this case there is dilepton emission 



3 



from current, medium, and constituent quarks, as well as from pion scattering with virtual 
bremsstrahlung. At T = T c nuclear matter transforms into an hadronic state, mainly of 
pions. For temperatures T c > T > Tf ~ 140 MeV, pion scattering with virtual quanta 
emission contributes. At the break-up temperature Tf we are left with only freely streaming 
particles. 



III. SOFT VIRTUAL PHOTON APPROXIMATION 

To begin the discussion, let us first recall some facts about bremsstrahlung in hadron- 
hadron collisions. If charged particles undergo acceleration, photons (real or virtual) may be 
radiated. Soft photon emission occurs only from the external lines [|16|| , as depicted in Fig |l|. 



Relative to this, radiation from the strongly interacting blob is negligible. It continues to be 
negligible as long as the energy carried by the photon is less than the inverse of the strong 
interaction collision time r ~ 1-2 fm/c; that is, if < 100-200 MeV. The advantage of the 
soft photon approximation is that the electromagnetic and strong interaction components 
of the matrix elements disentangle leaving separate multiplicative factors. This allows the 
invariant cross section for scattering and at the same time producing a soft real photon of 



four momentum g M to be written as [16 



<i 4 <r 7 a 
^° d 3 qdx 4tt 2 



pol 
X 



where da/dx is the strong interaction cross section for the reaction ab — > cd, e x is the 
polarization of the emitted photon, and 

J» = -Q a — - Q b — + Q c — + Q d — (3.2) 
p a ■ q Pb-q Pc-q Pd-q 

is the current. The Qs and ps represent the charges (in units of the proton charge) and 
the four-momenta of the particles, respectively. For soft virtual photons (dileptons) it is 
possible to continue from q 2 = to q 2 = M 2 , M being the invariant mass of the e + e~ pair. 
The result is [[[(J (q^ = p+ + p^t 



E E d 6 a e+e a 1 ^ d 3 a~< ^ 

d 3 p + d 3 p^ 2tt 2 q 2 d 3 q 

Eqs. (p.l|)-(p73|) can be combined to give the cross section for quasi-elastic ab — > cd scattering 
while at the same time producing an e + e~ pair. Four momentum conservation is strictly 
enforced at the photon-dilepton vertices by writing 

where t is the four-momentum transfer in the ab — * cd collision. We note that this is 



somewhat different from the expression used by Gale and Kapusta |L7] but this requires no 
further approximations whereas in theirs q$ was replaced with (E + + E_)J(E + + E_) 2 — M 2 
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in order to facilitate analytic integration. Such a replacement introduces a slightly different 
mass dependence into the result than does Eq. [3.4| . The expression from Ref. [17 is larger 
than Eq. |3.4| for small invariant masses and smaller for large invariant masses. A region of 
intermediate mass exists where the two are nearly the same. 

Separation of the full amplitude into strong and electromagnetic multiplicative factors is 
convenient. The strong interaction on-shell cross sections can be parametrized or modeled 
within, say, a quantum field theory. The electromagnetic factor can also be calculated 
exactly. It represents a coherent sum of all radiation fields involved, both initial and final 
charged states. Consider the completely general scattering process 



ab — > cd 



(3.5) 



where the masses m a = m c and m b = are arbitrary and the charges are subject only to 
overall charge conservation, namely, Q a + Qb = Qc + Qd- We calculate the exact squared 
modulus of the polarization dotted into the current to be 



J' 2 



ab^cd 



-(Ql + Ql + Ql + Ql) 



-2{Q a Q b + Q C Q 



s — m„ — mt 



y/[s- (m a - m b ) 2 ][s - (m a + m b ) 2 } 



x In 



m„ 



m b ) 2 + 



(m a + m b ) 



[m n 



+2(Q a Q c + Q b Q 



- m b ) 2 - 
2m a m b 
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x In 
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^Am a m b -t- \/^-t 



+2(Q a Q d + QbQc 



mt 



^[s + t- (m a + m b ) 2 ] [s + t - [m a 



m b y 



x In 



U + t- (m a - m b ) 2 + ^Js + t- (m a + m b ) 2 



>s + t- (m a - m b ) 2 - Js + t- (m a + m b ) 2 



(3.6) 



Details of its lengthy derivation can be found in Appendix A. The usefulness of this formula 
is in its generality. It is a benchmark with which to compare approximate expressions for 
|e • J\ 2 , which is our next task. If in this previous scattering process, the momentum transfer 
is small relative to the masses in the problem, specifically if \t\ < (m a + m b ) 2 , then a good 
approximation to the electromagnetic amplitude is (for details see Appendix B) 



I 71 ~- 
l C ' J lab^cd — 



2 1 



-t 



3 go \m a m b 



(Ql + Q 



;Q a Qbf(s) 



(3.7) 



where 



/(*) 



W- 

2^1 



2w_ 



(tU- — W+) J 2y 

4w- 
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where w- = s — (m a — nib) 2 and w + = s — (m a + nib) 2 . The interference function is shown in 
Fig. |2| for equal-mass scattering m a = nib = ni. Its behavior is an amusing consequence of 
multipole interference. In the limit \fsjm — > 2, f(s) — > —2/3 and as \fs/m — > oo, f(s) — > 0. 
Eqs. ( p.7|) and ( |3.8|) are extremely useful because the result is linear in t. It turns out that 
this approximation is quite good as long as particles are at least as massive as pions. For 
processes involving lighter particles, one must resort to Eq. ( |3.6| ). 

Whichever expression is used, either Eq. ( |3.6|) or Eq. (|3.7| ), the common l/qfi dependence 
allows analytic integration over the virtual photon and lepton momenta when inserted into 
Eq. Q3.4|) . The resulting differential cross section is 



da^cd _ g(f) L ( (V~s- K + m b )) + - K + m h )f - 

dM 2 4vr 2 M 2 | n \ M 



M 2 



\ {Vs- (m a + m b )) 



2 ( ' 

(3.9) 



where we have defined 



as 



/ dt^^ (q 2 \e ■ J\ 2 ab _ cd ) . (3.10) 



da, 

L_ 

-\(s,ml,ml)/s 



Sometimes called the kinematic triangle function, X(x, y, z) = x 2 — 2(y + z)x + (y — z) 2 |TJ 
Evaluation of the e + e~ cross section proceeding along these lines neglects the momentum 
of the virtual photon in the phase-space 5-function for the on-shell elastic cross section. 
Strictly speaking, one should correct for this by including the factor 

R 2 (s 2 ,m 2 a ,m 2 b) = X 1 / 2 (s 2 ,m 2 a ,m 2 b ) s 
R 2 (s,m 2 a ,ml) A 1 / 2 (s, m 2 a1 m 2 b ) s 2 

(with s 2 = s + M 2 — 2y/sqo) before integrating Eq. ( |3.4|) Hi~9"f . We will include this correction 
in the numerical part of this study. The independent particle approximation from kinetic 
theory allows the rate (number of reactions per unit time per unit volume) to be written as 



d^xdM 2 9ab J S J (27r)3 7 (2vr)3 e dM 2 Vrel ° V s ^ + P»> 

(3.12) 



where 



Vrel 



\f(Pa-Pb) 2 ~ ni\m 



E a E b 
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(3 is the inverse temperature and g a b = (2s a + l)(2s6 + 1)N^N^ is the spin and color degen- 
eracy. Integrating over both elastic scattering particles' momenta (for calculational details 
see Appendix C) yields 



dz 7Fi ^W-HTTW ( 3 - 13 ) 



d A xdM 2 16vr 4 J T 4 w rfM 2 

Zfn. i n 

where z m i n = (m a + rrib + M)/T, z = y/s/T, and JCi is a modified Bessel function. We 
are now in position to calculate the rate of electron-positron production through virtual 
bremsstrahlung from pion processes, quark-quark or quark-gluon processes, etc. 



IV. PAIR PRODUCTION FROM THERM ALIZED PIONS 



A truncated sum of the possible pion processes is used, that is, 



ab,cd \ / 



ab^cd 



+l<M*- ff ° 



da 
~dt 



\ € J Itt+tt- 



da" 
~dt 



(4.1) 



We refer to this as truncated because with it we are neglecting 7r + 7r + and 7r 7r processes. 
Such an omission is justified by the fact that there is significant cancellation due to de- 



structive interference in like-charge processes [pi. We have also neglected 7r + 7r — > 7r°7r° 



since a(7i + 7i — ► 7r 7r°) is rather small relative to the other pion processes | 21j| . Using the 
approximate electromagnetic amplitude, Eq. (|3.4|) becomes 



a gw(g) 

dM 2 ~ 6vr 2 M 2 



r(Vi,m.,M){^ + ^/(s)} 



(4.2) 



where a wn ~ cr-n+n + cr^-^o + cr w + w -, the momentum transfer weighted cross section is cx(s), 
and we have defined yet another function 



r( v ^,m 7r ,M) = In 



M 2 



M 



+ 



M 2 



\ {sfs - 2m 7 



,2 ' 



(4.3) 



If da/dt is a symmetric function of t and u, then JlO 

a(s) = 2aJs 



Am 2 



- 1 



(4.4) 



The elastic pion-pion cross section is parametrized in the following way: [|10 
a) for yfs < 0.6 GeV the chiral model expression is used, 
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2 11 

a e [[s) = - s 

v ; 3Fil6vr 



hmt 7m A 



1 ^ + 



s 2 



(4.5) 



with the pion decay constant F n = 0.098 GeV. 

b) at the collision energy near the rho mass, 0.6 < y/s < 1.5 GeV, the largest contribution 
to 7T7T the scattering amplitude is due to resonance formation. Therefore, 

/ \ _ g^nr (S - 4m 2 ) 2 
6 ' 1 J 487T S ( S -m2)2 +m 2 r 2' ^ 

where the coupling constant g p7T7T ~ 6, m p = 0.775 GeV, and T p = 0.155 GeV. 

c) for large collision energy y/s > 1.5 GeV, a e i becomes energy independent a e i ~ 5 mb. 
With a parametrization of the cross section we can combine Eqs. ( p,13| ) and ( (4.2|) to 

calculate the rate of e + e~ production from the sum of 7r + 7r°, 7r~7r° and 7r + 7r~ processes. 
Different regions of ^fs are investigated by varying the upper limit on the energy. The 
results at T = 200 MeV are presented in Fig. [5] along with a comparison of the rate using 
a quartic pion interaction Ci = \(tt ■ 7?) 2 /4 (A = 1.4 from tttc scattering lengths) p2| . Our 
results reproduce the rate from the quartic interaction in the low mass region M < 100 MeV 
only if a/s is restricted to the non-resonance region, say y/s < 0.48 GeV. Electron-positron 
pairs with invariant masses larger than ~ 100 MeV are attributed to yfs in the resonance 
region and beyond. 

Let us compare this production rate to a result which uses the exact electromagnetic 
factor. It is given by Eq. ( |3.6| ) with Q a , Qb, Q c , Qd = 0, ±1 and m a = = m n . Since the 
t dependence is no longer linear, knowledge of the energy dependent total elastic scattering 
cross section is not enough. We must have a parametrization or a model calculation of 
the differential cross section da e i/dt from which we integrate to obtain a e i(s). We perform 
a field theory calculation using a, p, and /(1270) meson-exchange to model the strong 
interaction. The quantum numbers for these particles are (I = 0, J = 0), (/ = 1, J = 1), 
and (/ = 0, J = 2) respectively. Low-energy nn scattering length analyses suggest that 



momentum-dependent scalar and tensor couplings to the chiral fields should be used |2o 
So the effective Lagrangian we use is 

A nt = g*<rdrf ■ c^tt + g p p " • (tt x «9 M 7?) + gff^n ■ cTtt. (4.7) 

Suppose we wish to calculate cr(7r + 7r~ — > 7r + 7r~). The differential cross section is proportional 
to the square of the matrix element describing the overlap of initial and final two-hadron 
states. There are six terms in this matrix element: t and s channel cr-exchange, p-exchange 
and /-exchange processes. The composite nature of the mesons necessitates modification at 
high momentum transfers. The net result is an effective suppression in that regime, so the 
vertices in the t-channel diagrams are given momentum-transfer damping monopole form 
factors 

h a {t) = (4.8) 
mi - t 



a 



where m a is either the cr, p, or the / mass depending upon which meson is exchanged. 

In this model finite resonance lifetimes are incorporated into the scalar, vector, and 
tensor boson propagators. For the / propagator we use |26[ 
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As previously mentioned, there are six terms comprising the matrix element. 



where 



Mi 
M 2 
M 3 

M 4 



M = Mi + Mi + M 3 + Ma + M 5 + Me, (4.10) 



-gihj(t){2<-ty 
t — ml + im a V a 

-gl{s-2mlY 

s — ml + im a T a 

-9 2 p h 2 p {t){s-u) 

t — m 2 n + im T, 



p 1 

9 2 p (u-t) 



: - m 2 p + im p T p 



(4.11) 



Notice that as written g a and gj are not dimensionless but g a mo and gftrif are. The full 
amplitude is symmetric under the interchange of the final states as required for bosons. The 
elastic differential cross section is 

^ = \Ml (412) 

dt 16vrs(s-4m2)' v ' ; 

It has qualitatively the right behavior: for energies just below the p resonance the angular 
distribution is forwardly peaked, for y/s near m p , the distribution becomes forward-backward 
symmetric, and finally above the resonance it is backward dominated. We can make a 
detailed comparison with data by integrating Eq. ( |4.12[ ) in order to get a e i(s). The scalar 
and tensor coupling constants and widths are adjusted so as to give resonant peaks which 
match experiment. Their numerical values along with the usual parameters for the p are 
g a m a = 1.85, m a = 0.525 GeV, T a = 0.100 GeV, g p = 6.15, m p = 0.775 GeV, T p = 0.155 
GeV, gftrif = 7.2, rrif = 1.274 GeV and Tf = 0.176 GeV. The results are shown in Fig. |] 



along with experimental data from Refs. [23] and [M. We obtain superb agreement with 



experiment for the energies shown. Other tttc processes (7r ± 7r° — ► 7r ± 7r ) proceed analogously. 
With all three cross sections, we can compute the rate of pair production using Eqs. ( p.9| ) 
and (|3. 10 ), which is exact within the soft photon approximation. In Fig || we show this rate 



along with the former (approximate) result of Eq. ( |4.2| ). The net result is a suppression from 
the approximate result by a mere factor ~ 2. Through this comparison we notice that e + e~ 
pairs with masses ~ 100 MeV are produced in tttt scattering at any kinematically allowed \fs 
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and are therefore dominated near the p-peak; whereas higher mass pairs (M ~ 300 MeV) 
are effectively produced through irn scattering with y/s ~ m p . When M ~ 300 MeV in 
Eqs. (|3.9| ) and (|3.13| ) with \fs near the kinematically determined minimum (just below the 
p-peak) the production rate is strongly suppressed by the logarithm and radical. Admittedly 
the soft photon approximation is less reliable for high masses but this suggests that it is 
imperative to have a reasonably accurate calculation of the cross section for y/s > m p . Such 
a calculation would be incomplete if it excluded the tensor /(1270) meson. 



V. QUARK PROCESSES 

Quark scattering processes introduce additional difficulties due to their fractional charge, 
their color, and their possibility for elastic scattering with their own force mediating bosons, 
the gluons. But we have developed the formalism to handle the general case. The same 
collection of diagrams from Fig. [l| is considered here with the charged external lines being u 
or d quarks (or antiquarks) and neutral external lines being gluons. In total we include six 
quark-quark (or antiquark) and four quark-gluon (or antiquark) diagrams. Other processes 
like qq — > gg and its reverse process turn out to give a rather small contribution to the e + e~ 
production rate (of order 8 % of our chosen processes) so we choose to neglect them here. 
It should also be mentioned that the approximate |e ■ J\ 2 cannot be applied to these. The 
electromagnetic amplitude for the partonic processes included is obtained from Eq. ( |3.6| ) by 
setting the masses to m q and by setting Q a and Q b to the appropriate quark electric charges. 
The strong interaction differential cross sections da qq /dt and da qg /dt are well known in the 



perturbative vacuum at the one-gluon-exchange level to be [£7 

dogb _ C ab 2^a 2 s 
dt t 2 ' 

where 



(5.1) 



For hot hadronic matter this is clearly inadequate. Let us consider the dominant t-channel 
gluon exchange amplitude for the diagram in Fig. |6|. A4 ab = agV^D^Y^, where T>^ u = g^ v jt 
is the gluon propagator and T a ^ b v are vertex functions for particles a, b. It is obvious that 
Ai ab is singular in the zero momentum transfer limit. But at finite temperatures static 
color-electric fields are shielded by quarks and gluons in the plasma, and color-magnetic 
fields are also shielded by non-perturbative effects . Color-electric shielding modifies the 
H = v = components of T>^ u in such a way that 

M% ^ a s (t)T a V 00 T b ~ ^rgl* (5.3) 

v 



(5.4) 



where the running coupling constant 

a E (t) a s (t) 



t t — ml 



E 
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and the color-electric mass m 2 E = 6na s T 2 . For a medium-modified magnetic scattering 
amplitude we obtain 



M ab ^ a s {t)Y1V^V) ~ _^L V ^ (5 5) 

t 



where the effective magnetic coupling etM{t)/t = a s (t) / (t — m 2 M ) and the color-magnetic 
mass m 2 M = ca 2 T 2 . Lattice gauge calculations give c ~ 20-30. In the small t limit the 
vertex functions reduce to ~ p" where is the momentum of particle a. Averaging over 
color and spin degrees of freedom gives the differential cross section |p8[ 



d(Tah „ 87T 



c, 



dt ~ ab ~sH 2 



aE(t)p a oP b -a M (t)p a -p b , (5.6) 



where s = (p a +pt) 2 - 

Unfortunately, da a b/dt depends on the quark-quark (gluon) orientation in the scattering 
through p a -pb- For arbitrary orientation of vectors p a , pb we can not obtain a closed expression 
for da ab /dt, but if vector p a is antiparallel to p&, da a b/dt reduces to 

da ab Tr a 2 (t)(2t-m 2 E -m 2 M ) 2 



dt a °2 (t-m 2 E ) 2 (t-m 2 



M) 



Note that this expression is exact for massless quark (antiquark) scattering. For T ~ 200- 
300 MeV, a s ~ 0.2-0.3, m E ~ 2.2T and m M ~ 1.5T. If we assume that the color-electric 
mass is roughly equal to the color-magnetic mass, m E ~ rriM, then we can restore the m q 
dependence in the differential cross sections 



da q9 2na 2 (t) 



dt s 2 (t-m 2 E ) 2K q 



<Y (5-8) 



da qq 4 27ra 2 (t) 



dt 9s 2 {t- m 2 E ] 



2K s-2mX (5.9) 



Using these expressions we did not observe sizeable changes in the dilepton production rates 
or final spectra. So instead we use Eq. Q5.71) henceforth. Note that in the limit the, tum — * 
(free space), we reproduce Eq. (O) for the perturbative vacuum. 



Putting together Eqs. (|3.4|), ( p.6| ), ( p,13| ) and ( f).7\ ), we have the rate of e + e~ production 
from all the quark(antiquark)-gluon or quark(antiquark)-quark processes. In Figs. [I] and 
H we compare usage of |e ■ J\ 2 approx to usage of |e ■ J\ 2 exact in calculating the rates of e + e~ 
production from the sum of qq{qq) and qg{qg) scattering with virtual bremsstrahlung at 
temperatures T = 200 and 300 MeV. Three quark curves are shown: current quarks, medium 
quarks, and constituent quarks with masses m q = 300 MeV. A remarkable feature emerges 
in the comparison between the correct and approximate rates from current quarks: there is 
a huge suppression. To perhaps better understand this effect we present in Fig. ^| the ratio 
R = \e ■ J\l xact /\e ■ J\ 2 approx for current, medium and constituent u and d quark scattering at 
y/s = 775 MeV. Owing to the small mass of the current quarks, the resulting ratio is R = 
10 -2 -10~ 3 . Another interesting feature for current quarks is that R — > 1 as cos# — > 1. Even 
at these (relativistic) energies, the electrodynamics in small-angle scattering is described 
very well by the approximate amplitude. However, these analyses invalidate Eq. (|3.7|) for 
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high momentum-transfer current quark scattering at these energies. Yet, for medium and 
constituent quarks the ratio R ~ 1, so for scattering of these particles the approximate 
expression, Eq. ( |3.7| ) is rather good. Fig. ^ also includes a curve for 7r + 7r~ scattering, where 
again the approximate electromagnetic factor is quite near the exact result. 

Since the quark mass dependence in the production rate is roughly 1/m 2 and the tem- 
perature dependence is roughly T 6 , an upper limit on the rate is obtained using m q = 5 
MeV at T = 300 MeV. Similarly, a lower limit on the rate is obtained using m q = 300 MeV 
at the phase transition temperature. Based on the results of Fig. [8] we conclude that at 
T = 200 MeV the (soft) virtual photon emissivities from competing medium quark versus 
pionic processes in the mass region M ~ 2m n differ by less than a factor 4. In a different 



energy regime Kapusta et al. [30] reached the same conclusion while comparing real (en- 
ergetic) photon emissivities of the QGP to that of hadronic matter. This equal-luminous 
property appears also in the annihilation channels qq — > e + e~ versus -k + tt~ — > p — > e + e~ 
but only after an (invariant mass) integrated total emissivity. In particular, the emissiv- 
ities dN e+e /d 4 x at T = 200 MeV for pions, medium quarks, and constituent quarks are 
approximately ~ 2.5, 1.2 and ~ 1.1 (x 10~ 9 GeV 4 ), respectively. 



VI. ZERO (TOTAL) MOMENTUM PAIRS 



As was pointed out by Kapusta |JT|, the production rate for soft dileptons could dif- 
fer by orders of magnitude from the naive prediction of quark-antiquark annihilation. The 
enhancement might be due to a modification in the quark dispersion relation in thermal- 
ized hadronic matter. Furthermore, he argued that the zero (total) momentum e + e~ mass 
spectrum is a direct measure of the quark dispersion relation and hence of chiral symme- 
try. A few possibilities for modifications are as follows. First, ignoring interactions with 
the matter would lead to a light-quark dispersion relation uj\ = k. Secondly, interactions 
with the plasma would result in a medium quark dispersion relation cu2 = (k 2 + T 2 ) 1 / 2 . Still 
another possibility is to use weak-coupling methods self-consistently [^], where one possible 



parametrization is lu 3 = (l/3 + k/T)k/(l + k/T). Having a modified relation between energy 
and momentum for the quarks leads to a modification in the annihilation rates for producing 
zero momentum lepton pairs. We can compare these rates with our soft photon approxima- 
tion only at a qualitative level since one calculation includes scattering with bremsstrahlung 
and the other includes quark-antiquark annihilation. Nevertheless, it is interesting to see 
how they relate. Another result with which we might compare more quantitatively is one 
which uses the resummation techniques of Braaten and Pisarski |33| . They were recently 
applied [jl2[ to a calculation of the zero-momentum soft dilepton production rate in a QGP. 



Here the differential rate for producing pairs is related to discontinuities in the photon self- 
energy. We refer the reader to Ref. [jl2[ for further details and merely state the results. The 



partial rate from the annihilation and decay of soft quarks and antiquarks exhibits structure: 
sharp peaks due to Van Hove singularities. In addition to the pole contributions, there are 
terms resulting from cuts in the effective quark propagator. The square of a coherent sum 
of such contributions appears in the final answer. The physical processes which give the 
pole-cut and cut-cut contributions involve soft and hard quark-gluon scattering processes, 
which we have considered in our bremsstrahlung calculation. At small invariant mass, the 
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pole-cut term grows like ~ 1/M 2 and the cut-cut term grows like ~ 1/M 4 ; they completely 
overwhelm the structure due to the pole-pole term. 

In order to compare our result with either calculation mentioned, we must restrict the 
lepton pair to have zero total momentum q = in the rest frame of the matter. With such 
restrictions in Eq. |3.4| , we cast the cross section into the form 

d A a e+e . a 2 a(s) 

where the weighted elastic cross section cr(s) is that which we have defined in Eq. BTTUL The 
rate dN/d A xdMd 3 q (q = 0) is then obtained by replacing da/dM 2 with d A a/dMd 3 q (q = 0) 
in Eq. [3 . 1 2| . In Fig. [10] we present the comparison of our (soft photon) approximate results, 
a parametrization of the results of Braaten, Pisarski and Yuan obtained using resummation 
techniques in QCD perturbation theory, along with the rate using the quark dispersion rela- 
tion co> 3 = (1/3 + k/T)k/(l + k/T) in the qq annihilation contribution. The other dispersion 
relations give even smaller rates. For dilepton masses M ~ 0.1 GeV, the production rate in 
the soft photon approximation is very near the result from QCD perturbation theory and 
much larger than the rate using the modified quark dispersion relation. Of course, the latter 
only includes annihilation whereas, the first two include scattering-so there is a qualitative 
difference. The low-mass dependence in the rates is ~ 1/M 4 both in the perturbation theory 
result as well as the soft photon approximation. This dependence is due to the processes 



of scattering hard quarks on hard gluons [12], which we have taken into account in the soft 



photon calculation. For masses 0.1 ~ M ~ 0.3 GeV, the ~ 1/M 2 dependence starts to 
become significant in the perturbation theory result. As previously mentioned, this depen- 
dence is due to the pole-cut contributions; or physically, with scattering of soft quarks on 
hard gluons. These processes can only be included in our scheme by retaining the next-to- 
leading term in the soft photon approximation. The matrix element M M in a soft photon 
approximation is expanded in powers of the photon momentum q as 

M> J ' = -^ + B l * + 0(q), (6.2) 

where both and B^ may calculated from knowledge of the physical, i.e. on-shell T- 
matrix |34j]. The term A^/q gives the 1/M 4 dependence in dN/d A xd 3 q (q = 0), while B^ is 
responsible for next-to- leading order terms. In particular only B^ ^ gives rise to 1/M 2 
dependence. The soft photon approximation with leading term only is quite close to the 
(resummed) perturbation theory result for masses M ~ 0.1 GeV. Most likely, the next-to- 
leading term (added coherently) would enhance the soft photon result for masses 0.1 ~ M ~ 
0.3 GeV by a factor of 2-6. However, such good agreement between our calculation and that 
of Ref. |12| does two things. It shows that scattering with virtual bremsstrahlung (rather 
than annihilation or Compton-like processes) accounts for most of the low-mass QGP-driven 
pairs and secondly, it lends support to the validity of the soft-photon approximation as 
applied to quark processes. 
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VII. TOTAL ELECTRON-POSITRON YIELD 



In the previous sections we calculated the rates of dilepton production by quark and pion 
scattering with virtual bremsstrahlung, i.e. the number of electron-positron pairs produced 
per unit four-volume d 4 x and with mass M while being locally thermalized at temperature 
T. In order to obtain the spectrum with which to make an experimental comparison, it is 
necessary to integrate over the time evolution of the nuclear system. Such an integration 
is possible only after modeling the evolution in some way. We shall make the standard 
assumptions |13|] about isentropic expansion and similarity flow in the central region (y = 0) 



of the heavy ion collision. Under such assumptions the initial temperature and time can be 
related to the hadron distribution in the final phase (for T > T c ) by 

T 3 r - ° 1 dNAA ^ (7 i) 

1 4a TtR\ dy ' 

where c ~ 3.6, a ~ 5.25, Ra is the nuclear radius (Ra = 7.1 fm for Pb), and Tj ~ 1 fm/c 
is the initial thermalization time. For LHC and RHIC energies the expected multiplicities 
dN AA ^ n I dy ~ 3000 correspond to initial temperature T ~ 300 MeV. As for the temperature 
of the phase transition T c , the lattice estimations provide a range of T c ~ 150-200 MeV. We 
assume that the phase transition is first-order, so at T c the mixed phase may be realized. 
For possible decay temperatures we use Tf = 140 and 100 MeV. The production in the pure 
quark and hadronic phases can be calculated || from 

dN 2 6 2 ?dT dm 



dydM A 1 1 J T 7 d 4 xdM 

and 

dN 2 6 2 2 } c dT dN* 

d# =3 ^ T ' r « r if«M' (7 - 3) 

while the mixed phase contribution from quark and pion processes are 



dydM 2 VfJ ^dHdM^ ' = ^ (7 ' 4 » 
and 

dN ^e d _ *Xk (Tl\\ Mr _ uJ^fT -T) (7 5) 

dydM ~ 2 \TJ T * r[r L) d*xdM [I ~ lc) - U ' 5j 

Here r denotes the ratio of the number of degrees of freedom of the QGP constituents to 
that of the hadron gas constituents (r ~ 12). 



In Fig. [XT] we present the resulting invariant mass spectra through virtual bremsstrahlung 
from quarks and pions using relevant temperatures Tj = 300, T c = 200 and Tf = 140 MeV. 
For this collection of temperatures the pion cooling phase contribution is the larger contrib- 
utor to its final spectrum, whereas the mixed phase dominates the resulting quark driven 
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spectrum. Note the difference between the results here compared to our initial estimates 
in Fig. 4 of Ref. JllJ. The suppression is due in part to the four-momentum conservation 
at the photon-dilepton vertices in Eq. ( |3.4| ) and to the phase-space correction shown in 
Eq. (|3.11| ). Also included in Fig. [TT] is an estimate of the background due to ir° and r\ Dalitz 
decay 7r° — > e + e~ + 7, 77 — > e + e~7. We calculate these backgrounds using formulae from 
Ref. This comparison suggests that 7r° Dalitz is the strongest source of pairs having 

invariant mass less than the pion mass and that rj Dalitz is the strongest source of e + e~ 
pairs with invariant mass larger that but less 2m v . For masses larger than 2m. w , tttc an- 
nihilation will of course becomes dominant. Superimposed on Fig. [TTJ is also the yield from 
QQ ~^ 7* ~~ > e + e~ annihilation in the Born approximation 0,|1(|. Throughout most of the 
range of dilepton invariant mass that we consider, the contribution from pion scattering with 
virtual bremsstrahlung and from it , rj Dalitz decays are two or three orders of magnitude 
above the qq annihilation spectrum. Perturbative corrections to this annihilation spectrum 
for finite p± are extremely important ||35|| . An advantage of our soft photon approach is that 
we could easily evaluate the bremsstrahlung contribution at finite q± or even impose some 
cuts in transverse momentum. 

Another quite reasonable possibility for the relevant temperatures of the colliding system 
is Tj = 250, T c = 150 and Tf = 100 MeV, corresponding to those used in Ref. [ID| in a 
recent study of low-mass dilepton production. We reproduce their tt^tt results using our 
|e ■ J\ 2 approx and using the same parametrization of the total titt cross section as did they. 
In Fig |12| we show our complete calculation using these temperatures. Our final results are 
rather different from Ref. [l(J since we include several improvements: we have other pion 
processes, we strictly conserve four-momentum at the photon-dilepton vertices, we use the 
exact electromagnetic factor and finally, we correct phase space to include the momentum 
of the virtual photon. It should also be mentioned that the pion results of Ref. |TU[ should 



not be taken seriously for invariant masses larger than ~ 2m n . They do not include any nn 
resonances in the cross section with masses greater than the p and most importantly they 
do not correct phase space which is essential for higher invariant masses. For comparison 
purposes, we show in Fig. 13 still a different collection of temperatures having the lowest T c 
which is compatible with lattice gauge theory data-Tj = 200, T c = 150, and Tf = 140 MeV. 
As one can see from Figs. |I^-[13|, the medium quark yield is always lower than the yields from 
7r°,77 Dalitz as well as pion scattering with bremsstrahlung. This being true, detection of 
quark degrees of freedom within this invariant mass region will be difficult in future heavy- 
ion experiments. Dalitz decays can (in principle) be distinquished from bremsstrahlung and 
removed from analyses. Distinguishing pairs produced through pion processes from those 
produced through medium quark processes will be difficult since total yields from both are 
nearly the same in all temperature scenarios we investigated. However, these experiments 
may be very important with regard to the study of pion dynamics in a thermalized system. 
For initial temperatures T; > 350 MeV {dN/dy ~ 4000 for LHC), pion bremsstrahlung is 
stronger than rj Dalitz in the mass region 100 ~ M ~ 300 MeV, but it is weaker than 7r° 
Dalitz for masses less than 100 MeV. From Eqs. ( [T.l| - [r.3D , we see the dilepton yield from 
thermalized pions is proportional to the square of the hadron multiplicity (dN/dy) 2 whereas, 
it is directly proportional to the multiplicity for Dalitz decays. Therefore, we expect that 
when the density of charged particles dN/dy ~ 4000, the yield of e + e~ pairs with 100 < M < 
300 MeV will be ~ (dN/dy) 2 , but with M ~ 100 MeV it will be ~ (dN/dy). 
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VIII. SUMMARY 



In a relativistic heavy-ion collision there is likely born a hot thermalized system. We 
argued that one of the best probes of such a system must surely be electromagnetic in 
its nature. In this study we have used low-mass electron-positron pairs as such a ve- 
hicle by investigating their origin, their abundance and behavior during these collisions. 
A soft photon approximation was used to compare hadronic to partonic scattering with 
virtual bremsstrahlung. In particular, radiation of virtual quanta during inr scattering 
with all isospin possibilities were discussed and compared with quark-quark(antiquark) and 
quark(antiquark)-gluon scattering accompanied by virtual radiation. We conclude that pi- 
ons and medium quarks are roughly equally luminous at T c , severely crippling our ability to 
separate one from the other. 

We compared our quark-driven production rate for zero momentum low-mass dileptons 
to that obtained using QCD resummation techniques. Our soft-photon-approximate result 
is very near the resummation rate for masses less than 100 MeV and a factor of 2-6 less 
than the resummation rate for masses larger than 100 (but less than 300 MeV). We expect 
that the next-to-leading term in the soft photon approximation (added coherently) would 
enhance our result for masses 0.1 ~ M ~ 0.3 GeV by 2-6. Such good agreement for masses 



less than 100 MeV between our soft virtual photon calculation and that of Ref. Jl2| using 



resummation techniques lends support to the validity of the soft photon approach to this 
problem. Using a (1+1) dimensional Bjorken hydrodynamic model to describe the expansion 
of the system, we calculated e + e~ total yields for different collections of temperatures T;,T C 
and Tf. For all reasonable sets of temperatures the constituent and the medium quark- 
driven total dilepton yields are lower than 7r° and rj Dalitz and lower than the pion-driven 
yield. We must conclude that it will be difficult to observe quark degrees of freedom (DOF) 
in this region of soft dileptons at LHC or RHIC. At the same time, it is important to keep 
in mind that our quark-driven soft-photon approximate results were below the pertubation 
theory results (shown in Fig. [10]) for invariant masses near 2m n by a factor ~ 4. This precious 
factor of 4 may be enough to allow reasonably clear observation of these DOF. Another useful 
observation is that for T- > 350 MeV {dN/dy ~ 4000 for LHC) the pion-driven production 
rate is larger than 77 Dalitz for masses 0.1 ~ M ~ 0.3 GeV, though still smaller than 7r° 
Dalitz for masses smaller than 0.1 GeV. Therefore, we predict that for dN/dy ~ 4000 the 
yield of dileptons at 0.1 ~ M ~ 0.3 GeV will be proportional to (dN/dy) 2 which reflects the 
collective effects in the thermalized pion gas. Whereas, the dilepton yield for M ~ 0.1 GeV 
is proportional to (dN/dy), i.e. tt° Dalitz dominated. Future heavy ion experiments (RHIC, 
LHC) will be useful for investigations of pion dynamics at high densities with regard to their 
soft dilepton mass spectra. 

In our discussion of pion processes we assumed that the tot differential cross section did 
not change from its free-function at finite temperatures. In particular, we used the same 
expression for da/dt and therefore o e \ near T = T c as we did for T = 0. It is clear that 
many-body effects introduce temperature dependences into the meson masses, widths and 
coupling constants. Therefore, the nir elastic cross section will likely be quite different near 
T c . In order to estimate the temperature dependence of da e \/dt and of a e i, we use as a first 



approximation the chiral perturbation theory modifications of F n and T p which are |36 
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F - (T) = F - (0) t 1 - im) • (8 - 1) 

and 

r p( T ) = i_TV4F2(0)- (8 ' 2) 

The zero temperature values are F n (0) ~ 0.098 GeV and r p (0) = 0.155 GeV. Since the 
pion decay constant changes, so does the coupling of the p to the pions because g pnw 
m 



/V2F n . The p-meson mass does not change appreciably with temperature ||37|| . We see 



from Eqs. ( | . 1|) and (|8.2| ) that the contribution of the non-resonance region (Vs ~ 0.6 GeV) 



in a e i will be enhanced with increasing T, whereas the contribution from the p-resonace will 
be suppressed. Basically, the p-peak gets widened and reduced with increasing T. It gets 
completely melted at T c . The spectrum dN e+e jdMdy obtained from Eq ( |3.13| ) by some 
portion of a space-time integration for bremsstrahlung during pion scattering is proportional 
to the area under the u e i{s) curve. Incorporating the temperature dependences in this fashion 
we observe little or no change in this area. This translates into the fact that we do not observe 
a large difference (less than a factor 2 enhancement) in the production rates or total yield 
with a temperature dependent cross section. Finally, we did not consider the change of 
T^o^g+g-^, or r\ ? _ >e + e - 7 because the lifetime of the 7r° and the r\ is such that only very few 
will decay during the mixed or hadronic phases. So the pion bremsstrahlung contribution to 
the dilepton spectra is not very sensitive to temperature effects, if temperature dependences 
are well accounted by Eqs. ( |8.1| ) and (^.2|). 

One topic somewhat ignored in our study is the transverse momentum dependence of the 
produced e + e~ pairs. Apart from the zero-momentum comparison we made, we presented 
all the results integrated over transverse momentum. Mostly because this study represents 
a first detailed analysis of quark and pion bremsstrahlung. There are at least two reasons 
we might benefit from transverse momentum analyses. Firstly, presentation of experimental 
data would mostly likely be g_i_-restricted due to detector (in) efficiencies and secondly there 
seems to be the possibility that Dalitz contributions might be suppressed more than would 
be bremsstrahlung for some restricted values of q± J39|. Finally, note that the soft photon 



approach that we have taken to this low-mass pair production problem can be equally well 
applied to the production of soft real photon production. 
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APPENDIX A: ELECTROMAGNETIC FACTOR |e J\ 2 

In this appendix we present the details for evaluating the exact electromagnetic factor 
for the soft-photon approximation. Recall the current in the scattering process ab — > cd: 

rryft nry^ 1 T)^ T)^ 

r = -Q a J ^- - Qb^- + Q c ^- + Qd^- (Al) 
Pa - q Pb- q p c - q Pd- q 



Summing over polarizations aided by current conservation leads to 3S 



^e^W = - J 2 . (A2) 

pol 

The particles' initial velocities are P a , Pb and final velocities are P c , Pd (c — 1). We move 
to the center-of-velocity frame of any charged participants where in the soft-photon limit 
\Pi\ = P for all i, Pa ■ Pb = Pc ■ Pd = —P 2 and P a ■ Pc = Pb ■ Pd = P 2 cos 9. Since p 2 = m 2 , we 
have 



|e- Jp 



QWa QW b Q 2 c m. 



E 2 a q 2 (l - Pa ■ n) 2 E 2 q 2 (l - p b ■ nf £ 2 g 2 (l - & • ff 



2 



QWd 2Q a Q b (l-P a -p b ) 2Q a Q c (l-p a -P c ) 



E 2 d q 2 (l - P d ■ nf q 2 (l - P a ■ n) (1 - P b ■ n) q 2 (l - P a ■ n)(l - & ■ n) 
2Q a Qd(l ~ Pa ■ Pd) 2QbQc(l ~ Pb ■ Pc) 



q 2 (l -p a -n)(l-p d - n) q 2 (l - p b ■ n)(l - /3 C ■ n) 
2QbQd(l - Pb ■ Pd) 2Q c Q d {l - p c ■ p d ) 



ql{l - p b ■ n)(l -Pd-n) q 2 {\ - P c ■ n)(l - p d ■ n) 



(A3) 



where n = q/\q\ is the unit vector for the photon. We are interested in an angular average 
over the emitted photon's direction 

^| e .J|». (A4) 

Carrying out the integrals for the first four terms in Eq. (|A]|) is elementary. The result for 
each is just the charge squared over q 2 . The remaining terms are readily evaluated using a 
two-parameter Feynman integral 



If dx 



ab J [ax + 6(1 — x)] 2 
We show, for example, integration of the fifth term in Eq. 

-2Q a Q b (l + p 2 



(A5) 



dVt q 1 



4vr ql(l-p a .fi){l-P b -n) 
-2QMI + P 2 ) JdxJ^ 



471 q 2 o 



1 - n ■ [P a x + P b (l - x] 
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r 1 
-2Q a Q b (l + I3 2 ) / dx—, 

i ql I 1 - n .x + 
i 



-2Q a Q b {l + (3 2 ) J dx^— 



P a x + p b (l - x) 

1 



q 2 (-A(3 2 x 2 + A(3 2 x + 1 - (3 2 ) 



-2Q a Q b (l + (3 2 ) 
-A(3 2 

-2Q a Q b {\ + (3 2 ) 



dx- 



A(3q { 



o % 
i 

dx 

o 



1 1 



1 , I-/? 2 



+ 



X 2 + V 4 + 4/3 2 



1 I I-/? 2 



(A6) 



4/3 2 



The integration is then trivial; yielding, 



. 2QaQt(1+ ^_l_ ln (I±|* 



(A7) 



Evaluation of all other terms in Eq. (|A3|) follows the same prescription. The net result of 
such an exercise is 

\e ■ A 2 = \l - {Ql + Ql + Ql + Ql) - 2{Q a Q b + Q c Q d ) i±£ In ( ' ' ^ 



+2{Q a Qc + QbQo 



2(3 



7(1 - (3 2 cost 



(3yJ(l - cos0)[/?V(l - cos 9) + 2] 



x In 



/3 2 7 2 (1 - cos#) + 2 + v//3 2 7 2 (l - cos 



+2(g a g d + g b g 



^//3 2 7 2 (l - cos 0) + 2 - 0? 2 7 2 (1 - cos 1 
7(1 + /5 2 cos#) 



/3^/(i + cos^)[/?V(l + cos#) + 2] 



/32 7 2(l + cos^) + 2 + x //5 2 7 2 (l + cosi 
//3 2 7 2 (1 + cosfl) + 2 - v //3 2 7 2 (l + cost 

The velocity and angle can be replaced with invariants through 



(A8) 



s - (m a + m b f 
\s- (m a - m b y 



(A9) 



and 



cos 9 = 1 + 



2f 



(m a + m b y 



(A10) 



In future calculations this general expression will be very useful, but for our purposes here 
it is enough to use equally massed particles, m. Then Eq. (|A9|) reduces to 
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and Eq. ( |A10p reduces to 



/ s — Am 2 

P = \ — . ( AU ) 



2t 

cos# = l + — -. (A12) 

s — Am z 



Carrying out these substitutions in order to express the result in terms of invariants s and 
t we arrive at 

\fJ\ 2 =\{-{Ql + Ql + Ql + Q 2 d ) 

*>fr> n . n n \ s ~ 2m2 , ( V$ + Vs-Am 2 \ 
-2(Q a Q b + Q c Q d )^= In — 

s(s-4m 2 ) Vv^- y/s-4m 2 J 



\s + t-Am 2 ){s + t) Vs + t-Am 2 , 



2m 2 -t ( J Am 2 -t + \f = i\ , 
+2 Q a Q c + Q fc Qd . In V = = V , = , A13 

ft(t - Am 2 ) V V 4m 2 - t - y/^t' 



where y/s > 2m and — (s — 4m 2 ) < t < 0. 



APPENDIX B: APPROXIMATE ELECTROMAGNETIC FACTOR 

It is perhaps more convenient to have an approximate, yet reasonably reliable expression 
for |e ■ J| 2 . This appendix is devoted to such a calculation. Let = (E a ,p a ) be the four- 
vector for particle a, and similarly for particles b, c and d. Let = qo(l,n) again be the 
photon's space-time description and choose the radiation gauge so that = (0, e). Then 
the quantity |e • J\ 2 can be written 



\e-J\ 2 = - 2 



QcPc 



QaPa 



+ 



QdP d 



QbPb 



(1-H-P B ) (l-n-Pa) (l-n-(3 d ) (l-n-/3 6 ) 



(Bl) 



Let us work in the center-of-velocity frame of any charged participants and write j3 a = /3, 
Pc = P + A/3, Pb = —P and P d = — (/3 + A/3). We also assume that Q a = Q c , Qb = Qd- Then 
if the scattering is such that A/3 is small relative to /3, a good approximation to Eq. (|B1|) is 



4 



'Q a (Ap + nx (px A/3)) Q b (-A/3 + n x x Ap] 



[1-n-p) 2 



(i + n-py 



(B2) 



where /3 is the initial (or average) velocity. Since we consider small angle deflection, A/3 is 
approximately perpendicular to the incident direction. Next we explicitly sum over polar- 
izations by defining the unit vectors eii and e± to be the polarization vectors parallel and 
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perpendicular to the plane containing (3 and n. Averaging over the azimuth of the scattered 
particles, the two polarization summands become 



2= \APl rdQ(_ 
2ql J An \{l 



Qa 



+ 



and 



\e-J\ 



2qi 



pcosO) (l + ^costf), 



<m f Q a ([3-cos6) Q b ({3 + cos6) ' 
An~ V(l -f3cos$) 2 + (1 + (3cos6) 2 



(B3) 



(B4) 



Integration of these expressions is elementary but quite lengthy. We perform the necessary 
work and add the contributions from both polarizations to get 



\e-J\ 



|A/3|V f Ql , Ql , (Ql + Ql) , 



approx 



<i 



QaQb 



2 7 2 



[66 

1 1 (3 
- H h — 

/5 2/? 3 2 



1 QaQb 

(3 2 2 



In 



1-/3 



(B5) 

For arbitrary masses m a = m c and m^, = rrid, 7 2 |A/3| 2 = —tj{m a m^) and /3 2 = {s — (m a + 
m fc) 2 }/{ s — ( m « — n^b) 2 }- It is enough for our purposes here to use equally massed particles 
m a = m b = m, in which case 7 2 |A/3| 2 = —t/m 2 and f3 2 = (s — Am 2 )/s. Restoring the 
invariants s and t, we get the final (approximate) result of 



\e-J\ 



approx 



2 1 

3 q 2 \m 



-t 



(Ql + Ql) - 7;QaQ b f(s) 



(B6) 



where 



Am 2 



2(s - 4m 2 ) 



m 
s 



2s 



+ 



V s — Am 2 \ \Js — Am 2 



V s — Am 2 1 j / y/s + Vs — 4m 2 
J \ - y/s - 4:m 2 / 



(B7) 



Strictly speaking, this expression is valid only for \t\ < Am 2 



APPENDIX C: RATE CALCULATION FROM GENERAL SCATTERING 

PROCESS 

The starting point for the rate calculation (from kinetic theory) is 

= - / - / / (- - <* + »>') . 



(CI) 
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where g ab = (2s a + l)(2s 6 + 1)N^N^ is the (spin and color) degeneracy. Trivial angular 
integrations can be done 

= 32^J dS J dPa J dPb E a E b dM* X ^ m - m ^ ^ 

where 

A(s, m 2 a , m\) = s (s - 2{m 2 a + m 6 2 ]) + (m 2 a - m\f (C3) 

and 

X = [4£ a £ fe (s - m 2 - m 2 ) - (s - m 2 - m 2 h f - Am 2 a E 2 - Am 2 b E 2 a + 4m 2 m 2 ] . (C4) 
Using the substitutions 

u = E a + E b and v = E a - E b , (C5) 

the differential element changes, due to the fact that the Jacobian is different from 1, as 
follows 

dE a dE b -> (C6) 



Completing the substitution, we get 

JdsJdvJ due-^^^X^^mlml), (C7) 



d 4 xdM 2 64tt 4 
where 

u± = ± ^XX)^ v /^_ (cg) 

After integrating over -u, the remaining v integration can only be written in terms of a 
modified Bessel function. The rate is 

d^xdM*- 32tt 4 J d8 ^ m ^ m b) pV ~ s dM 2 ■ 
Introducing the dimensionless variable z = y/s/T, we arrive at the final expression 

d^dM 2 = 16^ J dZ T 4 " ~dM 2 ~ (C10) 

where z min = (m a + m b + M)/T. 
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FIGURES 

FIG. 1. Lepton pair production through virtual bremsstrahlung in scattering of particles 
a b — » c d. The particles might be charged or uncharged pions, quarks or antiquarks or gluons. 
The shaded region indicates a strong interaction. 



FIG. 2. The interference function f(s) as it changes with the dimensionless variable yfs/m. 

FIG. 3. (Approximate) e + e~ production rates through virtual bremsstrahlung from pions as 
compared with quartic pion interaction predictions. The dashed lines are our estimations with 
varying values for the upper limit on the energy and the solid line is from the quartic interaction. 



FIG. 4. The total elastic 7r + 7r cross section as compared with experimental data from Refs. [23 



and |4|. 



FIG. 5. e + e~~ production rates from pions using the exact electromagnetic factor and the 
model-calculated differential cross section (solid curve) as compared to the rate using the ap- 
proximate electromagmetic expression and the assumed-symmetric parametrization of the total 
cross section (dashed curve). 

FIG. 6. i-channel gluon exchange diagram. 

FIG. 7. Quark-driven rates at temperatures T = 200 and 300 MeV along with the pion driven 
rate at T=200 MeV using the approximate electromagnetic factor applicable to each process. The 
long-dashed curves are for current quarks, the dot-dashed curves are medium quark results and 
the short-dashed curve (also at T = 200 MeV) is for constituent quarks. 

FIG. 8. Quark-driven rates superimposed on the pion-driven rate but this time using the correct 
electromagnetic factor. The labelling scheme is the same the previous figure. 

FIG. 9. The ratio R = |e • J\exact/\ e ' J\approx as ^ depends on the center-of-mass scattering 
angle cos 9. 

FIG. 10. Production rates for zero (total) momentum dileptons: The solid curve is a 
parametrization of the results from Ref. |i"2|| , the dashed curve is our soft-photon approximate 
results and the dotted curve is the result one gets by modifying the quark dispersion relation. 

FIG. 11. The total e + e~ yield from competing sources. The solid curve is for the pion-driven 
yield, the short-dashed curves are Dalitz decays, the long-dashed curve is the Born approximate 
qq annihilation result and the dot-dashed curves are I) current quarks, II) medium quarks and III) 
constituent quarks. 
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FIG. 12. The total e + e yield with different relevant temperatures. 
FIG. 13. The total e + e~ yield with still a different set of relevant temperatures. 
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